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Abstract 

An analytical model, recently developed to study the electromagnetic dissociation 
(EMD) energy spectra of two-neutron halo nuclei, is applied to ^^Li and ^^Be. 
We find that a reliable set of experimental data on EMD could help to resolve 
the problem concerning the structure of the ^^Li ground state. For ^^Be we find a 
mutual inconsistency between the existing experimental data on the binding energy, 
radius and EMD energy spectrum. We also conclude that the structure of ^^Be is 
essentially more complicated than the structure of other Borromean two-neutron 
halo nuclei. 
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1 Introduction 



The reaction of Coulomb excitation in collisions with heavy targets is an 
important tool to study the structure of unstable nuclei. For nuclei having 
only one bound state the excitation means breakup and the process is called 
electromagnetic dissociation (EMD). In the present paper El EMD of two- 
neutron halo nuclei is considered. Microscopic calculations of the process at 
the three-cluster level were performed for ^He [1,2], ^^Li [2,3], and ^^Be [4,5]. 
Such calculations are hindered by incomplete knowledge of the cluster dy- 
namics. Quantum Monte Carlo A-nucleon calculations [6] are hindered by 
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the complexity of the A-body problem. Therefore it seems to be of interest 
to describe the process semi-phenomenologically aiming both to understand 
its main features and to guide experimental studies. Such an approach was 
developed in our recent paper [7]. Three-cluster bound-state wave functions 
(WFs) of two-neutron halo nuclei were constructed which behave correctly at 
large intercluster distances, reproduce the nuclear sizes, and incorporate the 
main features of the underlying three-body structure. No-interaction three- 
body WFs were used to describe the breakup final states. The model allows 
an analytic calculation of the strength functions and can serve as a tool to 
predict the Coulomb disintegration spectra of a variety of two-neutron halo 
nuclei. It may be used when precise knowledge concerning the structure of 
a system under study is lacking. As a test case the relatively well-studied 
nucleus ^He was considered in [7]. We found that the large-distance asymp- 
totics of the ground-state WF determines the shape and the position of the 
maximum of the strength function, while the size of the ground state governs 
the asymptotic constant and thus the magnitude of the strength function. As 
to the final-state interaction (FSI) effects, one can note that the total El 
strength will not change when replacing a complete set of true final states 
with a complete set of no-interaction final states. Therefore the inclusion of 
FSI will result only in a redistribution of the strength leading to a somewhat 
higher strength at low energies. But even without including FSI a very good 
agreement concerning the shape and peak position of the El strength function 
and a reasonable agreement concerning its magnitude was found in the ®He 
case [7]. 

In the present paper we apply the model to study the Borromean halo nu- 
clei ^^Li and ^'^Be. For both nuclei there are large ambiguities concerning the 
ground-state structure. Furthermore, there are large uncertainties in the ex- 
perimental data regarding EMD reactions and, for ^^Be, also concerning basic 
ground-state quantities such as binding energy and size. Therefore we believe 
that an analytical model with large freedom in trying different assumptions 
can be very useful. 

In Sec. 2 the model is outhned, more details can be found in Ref. [7]. In Sec. 3 
the results of our analysis of the ^^Li and ^^Be EMD data are presented. Sec. 4 
contains the conclusions. 



2 Analytical El strength functions of two— neutron halo nuclei 

We consider the electromagnetic excitation of a system with only one bound 
state and consequently with all possible final states belonging to the contin- 
uum. In the framework of first-order perturbation theory the energy spectrum 
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for El Coulomb excitation can be written as 

d(7(El) _ iVEi(^*) 167r3dfi(El) 
dE ~ he 9 dE 



Here E* is the excitation energy, E* — Eq + E, where £'0 is the binding 

energy, and E is the continuum energy, Nei{E*) is the spectrum of virtual 
photons [8,9] and dB{El)/dE is the dipole strength function 

^^ = ;t7Vt^ ^ [ dTf\{f\M{El,n)\t;JM)fS{Ef-E), 

(2) 

where dr/ is the phase-space element for final states, M.{E1) is the dipole 
transition operator, \i) is the initial state, and |/) are the final states in the 
center-of-mass (CM) system. The notation / drj in Eq. (2) implies summa- 
tion over discrete quantum numbers in addition to integration. We will study 
halo nuclei which are known to exhibit a large degree of clusterization. Since 
we are interested in low-energy excitations, involving mainly relative motion 
between the clusters, we will use the corresponding A^-cluster El operator 

M{El) = 1/ — E ^Ziin - Rcra), (3) 

V 47r 

— * 

where fl are the cluster positions, and i?cm is the position of the CM of the 
system. 

We adopt the three-body model to describe two-neutron halo nuclei. The clus- 
ter part of the bound-state WF, in the CM system, is written as an expansion 
over hyperspherical harmonics (HH), see e.g. [10], 



^{x,y)=p-'/' E XkIs{p)[^kI{^^)^Os]j^. (4) 
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Here {x, y] is the set of normalized Jacobi vectors where x corresponds to the 
vector joining the two valence neutrons, and {p^VL^} are the hyperspherical 
coordinates in the {x, y] space. The functions T^'f^ are HH, and Osms {S — 
0, 1) are the spin functions of the two valence neutrons. Since the HH expansion 
converges rapidly we retain only one, or a few, terms. Furthermore, since the 
hyperradial functions, Xkls{p)^ should behave rather similarly at large p, 
which is the region of interest, we use the same hyperradial function for all 
terms retained in the expansion. These approximations lead to the following 
normalized WF for initial bound states having J'" = 0"*" 



+O40r4oo (^5) ^00 + Ct21 ^21 



(5) 
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The coefficients aoo, 0-20; 0,21 and a^Q entering here are the amphtudes of various 
HH, Oqq + + ^21 + '^40 ~ addition to the WF components used 

in [7], the K = 4 contribution is introduced here and in Eq. (9) below.) We 
consider hyperradial functions of the following two forms. The first is a simple 
normalized exponential 



'{p) = V2Koexp{-Kop), (6) 

where the single free parameter is fitted to the binding energy via Eq = 
{hK^Y / {2m) ^ and m is the nucleon mass. This function reproduces the true 
asymptotic behaviour of Borromean three-body systems with two neutrons as 
constituents, see e.g. [11]. However, the three-body size (p^), and consequently 
the total size of the system, is underestimated with such a model. Therefore, 
we introduce a second hyperradial model function having two free parameters 

X^^Hp) = c [exp(-/«op) - exp(-ft;ip)] , 
where c 



The parameters Kq and Ki are fixed using experimental values of binding en- 
ergy and rms radius, see [7]. The condition that ni > k.q ensures that the 
second term decays faster than the first, and thus the correct large-p asymp- 
totics is preserved. 

We use a complete set of no-interaction final states that include just coordinate- 
space HH 



^Kim^Os]^^. (8) 



Alike the plane waves, these states are solutions to the free-space six-dimensional 
Schrodinger equation. The continuum energy is related to k via E = {hnY / (2m) . 

Using these model WFs an analytical calculation of the transition matrix 
elements (MEs) entering Eq. (2) is possible. When the model WF of Eqs. (5, 
7) is adopted for the bound state, one obtains the following final expression 
for the El strength function: 



,%m + E){Ej + E)f/' 
X [aiFi{yi)Fi{yj) + a2F2{yi)F2{yj) + a^F^{yi)F^{yj)\ , (9) 

Here £'0,1 = {fiKQ^if / {2m), where Kq, «i and c are defined in Eq. (7). The 
constant D is 

2l^2my {A-2)A' 
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and the constants cti, 0:2 and 0:3 are 




9009v^\ V 2 2 4 2 4 \ 

7^ ( ^20 + «21 + n«40 + Q «20«40 , 



(4aoo + 020 + «2i + 4000^20) , 
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The coefficients ooo, 020, 021 and 040 are defined in Eq. (5). The quantities 7/0,1 
are defined as yi — E/(E + Ei), and, finally, 



where F{a, (3, 7, z) is the standard hypergeometrical function. The correspond- 
ing three terms in Eq. (9) represent contributions from final states with K — 1, 
K = 3 and K = 5 respectively. When the hyperradial function (6) is used the 
result is obviously obtained from Eq. (9) by retaining only the term with 
i = j = and replacing with 2k:o. 



3 Results 



Despite the great number of experimental and theoretical studies that have 
been performed, the ^^Li nucleus still remains a puzzle. The role of possible 
s-wave intruder states [12-16] is of much theoretical interest. The experimen- 
tal evidence for virtual s-states in the binary subsystem ^°Li is diverse and 
sometimes confiicting [17-24]. Three-body microscopic calculations of ^^Li are 
uncertain due to this fact, in particular. Note that whereas the most recent 
experiments establish that the ground state of ^''Li is a si/2 state [24], this it- 
self does not mean that the configuration is a predominant one in ^^Li. 
The weights of the (si/2)^ and (pi/2)^ configurations depend on the relative 
strength of s- and p-wave potentials determining the gap between the Si/2 
and pi/2 states in ^°Li. The shape of the ^°Li momentum distribution, mea- 
sured in a ^^Li fragmentation reaction [25], supports a (lsi/2)^ contribution 
of (45 ± 10)% in the ^^Li ground-state WF. Further evidence for the mixing 
of different parity states is clearly seen in the asymmetric angular correla- 
tions [25]. 




3.1 
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Li 
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Clearly, the present-time understanding of the structure of ^^Li is substantially 
poorer than that in the ^He case, in which case more definite information on 
the interaction between the valence neutrons and the core is available. This 
makes ^^Li an interesting nucleus to study with our analytical model. In this 
way we will be able to test different assumptions on the ground-state structure. 
In the present work we have used four different ground-state WFs of ^^Li, 
which are summarized in Table 1. Here we restricted ourselves to HH with the 
lowest possible K values, K = Q and K = 2. These configurations correspond, 
respectively, to the (si/2)^ and configurations of the core-centered shell 

model for valence nucleons. 

To explain the latter fact we need to consider a relationship between HH and 
the configurations of the core-centered shell model. One can pass from the Ja- 
cobi vectors {x, y} defined in connection with Eq. (4) to another pair {xi, yi} 
of Jacobi vectors, such that Xi is proportional to the distance between one 
of the valence nucleons and the core, and yi is proportional to the distance 
between the other valence nucleon and the CM of the pair consisting of the 
former valence nucleon and the core. The distance from the CM of this pair 
to the core is times smaller then that to the valence nucleon belonging 
to the pair. In the framework of the approximation where one neglects the 
former distance as compared to the latter one, i.e. locates the CM of the pair 
on the core, the Jacobi orbital momenta associated with the vectors {xi} and 
{yi} coincide with orbital momenta of the two valence nucleons in the core- 
centered shell model. In order to find HH to be retained in the model WF 
we proceed as follows. Let v''^lj^{x,y) denote the HH entering Eq. (4). In 
addition, we define the HH '^^KLMi^i^Vi) where l^i and ly^ are Jacobi orbital 
momenta associated with xi and yi. We consider first the expansion over the 
latter HH and we retain the terms in this expansion for which the orbital 
momenta L, l^^., and ly^ are the same as for predominant configurations in the 
core-centered shell model. Note, however, that the corresponding HH do not 
coincide with the lowest shell-model configurations but also inchidc higher 
shell-model states with radial excitations. The final formula (5) with the coef- 
ficients listed in Table 1 is obtained expressing the retained HH '^^KLMi^iiVi) 
as linear combinations of the HH T ^l}^^{x^y). The quantum numbers K and 
L are conserved at this transformation of HH. 

The first two of our WFs are 100% K — Q states. In the case oi K — states 
the orbital momenta of relative motion in each pair of the three particles 
are zeros. Thus the K = Q state corresponds to the (si/2)^ configuration of 
the core-centered shell model. The wave function \l'i(^^Li) contains the one- 
parameter hyperradial function (6) while ^2(^^Li) contains the two-parameter 
function (7). 

The two other WFs, *3(^^Li) and ^4(^^Li), include HH with = and 
K — 2. They correspond to 50% mixtures of (si/2)^ and (pi/2)^ configurations. 
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To obtain the K = 2 HH and their weights hsted in Tabic 1 wc expand the 
J = neutron shell-model state over such states in LS coupling. 
In the framework of our infinite core-mass approximation the latter states 
correspond to HH, T^'^^, with K — 2 and li — I2 = I, defined with respect 
to the above mentioned Jacobi vectors {xi, yi}. Two such HH arise, one with 
L = 1, and the other with L = 0. Finally, when the above mentioned rotation 
to the {x, y} system is performed the first HH (with L = 1) transforms into 
a single HH with K = 2, L = = ly = 1 and the second one (with L = 0) 
transforms into a single HH with K — 2, L — = ly ~ 0. 

The WF ^4(^^Li) differs from ^3(^^Li) only in that the relative signs of the 
amplitudes pertaining to the (pi/2)^ configuration and the (si/2)^ configuration 
are different. 

The results of our calculation are presented in Fig. 1 along with available 

experimental data. The El strength functions arc plotted in the figure. In 
Fig. 1(a) the results are presented in arbitrary units and in Fig. 1(b) they 
are plotted in absolute scale. There exist three published sets of experimental 
data on EMD which show large ambiguities. The experiments were performed 
at different energies, Zinser et al. [26] used a high energy beam (280 MeV/A) 
while Sackett et al. [27] and Shimoura et al. [28] had much lower energy, 
28 and 43 MeV/A, respectively. The validity of extracting the El strength 
function from the low-energy experiments [27,28] is discussed in [27]. Sack- 
ett et al. report their data only up to ii^ = 1.45 MeV above particle-decay 
threshold, and Shimoura et al. list their strength function only in arbitrary 
units due to the lack of a forward-angle detector. In the Zinser et al. experi- 
ment, the strength function was decomposed into two Gaussian components. 
The two-peak structure observed by Zinser et al. was not seen in the other 
two experiments. It is clear that the data sets are contradictory. The shapes 
are different, and the peak amplitude is much higher in the Sackett et al. data 
than in the Zinser et al. data. However, in all data sets the position of the 
maximum is approximately the same. 

Our results may be commented as follows. All our ground-state WFs give the 
peak position at about 0.5 — 0.7 MeV which is very close to all experimental 
findings. This value is expected. Indeed, for not too high energy the predomi- 
nant contribution to the strength function comes from large distances, and in 
our model it is provided by the first, longer range exponential in Eq. (7), see 
also [7]. Because of this the strength function is mainly a function of E/ Eq. 
Moreover, in the peak region the dominating final state has K — 1 which cor- 
responds to the contribution of the single term with Fi in Eq. (9). Due to these 
facts the positions of the peak are close to each other for all our ground-state 
WFs and proved to be about 2Eq, see also [29]. However, terms with larger K 
in (9) change the behaviour at higher energies. There is no possibility to ex- 
plain the two-peak structure of Zinser et al. in our model. The one-parameter 
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WF, \E'i(^^Li), leads to a very low strength which is expected since it underes- 
timates the size of the system, see Table 1. The results obtained with our three 
other model ground-state WFs of ^^Li are to be compared with experiment if 
the absolute values of the strength function are considered. 

We see that the ^2(^^Li) WF, which corresponds to s-motion of the va- 
lence neutrons with respect to the core and also to relative s-motion be- 
tween them (or a configuration in core-centered, shell-model coordi- 
nates), compares rather well with the Sackett et al. data. However, for energies 
E > 0.5 MeV the error bars in [27] are large. On the contrary, the ^3(-^^Li) 
WF, where we have equal weights of (si/2)^ and (pi/2)^ configurations in ^^Li, 
leads to an excellent agreement with the shape of the Shimoura et al. data 
and a relatively good agreement with the absolute values of the Zinser et al. 
data. Finally, the ^4(^^11) WF leads to results which are significantly different 
from those obtained with ^3(-^^Li). This may be related to the fact that the 
corresponding difference in the relative sign of the (>Si/2)^ and (pi/2)^ compo- 
nents of the WF leads to a change of the internal geometry of the system, and 
the Tc value in particular, as it is seen from Table 1. Thus there exists a sen- 
sitivity to the relative sign of the two components in the WF and not only to 
their weights. In conclusion we state that this kind of analysis, when applied 
to a reliable set of experimental EMD data, would give valuable information 
concerning the relative role of the (si/2)^ and configurations in the ^^Li 
ground state. 



3.2 ^^Be 



The existing experimental data on ^^Be are scarce. Even such a fundamental 
quantity as the binding energy is characterized by large uncertanties. Com- 
bining the two published measurements gives the tabulated two-neutron sep- 
aration energy of 1.34 ± 0.11 MeV [30,31]. The most recent value of the rms 
radius is 3. 10 ±0.15 fm [32]. Experimental and theoretical efforts are presently 
focused on understanding the ground-state structure. Assuming ^^Be to be a 
nucleus where the neutron p-shell is closed one gets 0^5/2 as the next orbital 
for the valence neutrons. Several experiments have indeed confirmed a state 
around 2 MeV in ^"^Be which is assigned having spin-parity I'" = 5/2"*" [33-36]. 
However, in order to reproduce the large size and narrow momentum distri- 
butions, a large s-wave component is needed in the ground state of ^"^Be [4,5]. 
The s-wave potential fitted in [4,5] to reproduce the ^^Be binding energy, with 
the known position of the 0^5/2 state, leads to an s-wave bound state in ^^Be 
which is clearly not observed. This problem may be connected with the fact 
that the ^^Be nucleus has no closed neutron p-shell (see [37,38]) and this ef- 
fect has not been taken into account in [5]. Recent experiments still indicate 
the presence of a low-lying virtual s-wave state in -"^^Be [36]. Furthermore, 
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studies of the isobaric analog state of ^^Be give an s-wave spectroscopic factor 
of 45 ± 20 % in the -'^'^Be ground state [39] which is supported by conclusions 
made in other experiments [32,40]. The K = component of our model WF 
provides s-motion between all the particles, and we used the K — A compo- 
nent to include some d-motion. However, since the K — A component leads 
to a very broad distribution of strength the shape will mainly be determined 
by the K — component. 

The EMD energy spectrum for ^^Be has recently been measured by Labiche 
et al. [40]. This experiment was performed at 35 MeV/A. In their analysis the 
energy spectrum was modeled using the Breit-Wigner shape with parameters 
Ebw = 1-8 ± 0.1 MeV and Tbw = 0.8 ± 0.4 MeV, see Fig. 2 in [40]. We prefer 
to compare theory and experiment using strength functions which facilitates 
comparison with EMD cross sections for other nuclei and at other beam en- 
ergies as well. We therefore extracted the dipole strength function from the 
EMD spectrum of Ref. [40] dividing out the virtual photon spectrum entering 
Eq. (1). This results in a very narrow distribution peaked just below 2 MeV as 
can be seen from Fig. 2. Note that the shape of this strength function is very 
different from those of strength functions of other Borromean halo nuclei. It 
should also be mentioned that in the experiment [40] only ^^Be and neutrons 
were detected. Therefore, we should consider the strength function up to the 
^^Bc threshold at 3.2 MeV because above this energy other reaction channels 
are open. 

In a first step of the analysis we use the position of the peak to check the 
binding energy. If the ground-state WF has an appreciable A' = weight then, 
as mentioned above, the peak should be positioned at approximately 2Eq. As 
can be seen from Fig 2(a) the experimental spectrum would correspond to a 
binding energy of Eq fa 0.9 MeV, considerably lower than the tabulated value 
of 1.34 MeV. In fact, we have one more indication that the tabulated binding 
energy is too large. Namely, the binding energy seems to be incompatible with 
the radius of ^"^Be. In Table 2 we compare the hyperradii 



of ^He, ^^Li and ^^Be which represent the three-body sizes of the systems. 
They are calculated from the rms radii of these nuclei and their three-body 
constituents. For halo nuclei the hyperradii are mainly determined by the 
asymptotic behavior of the WFs. Therefore their values should exhibit a cor- 
relation with binding energies. Due to this reason one would expect from the 
systematics of ^He and ^-"^Li that either the binding energy of -"^^Be should be 
smaller at given size or the size should be smaller at given energy. 

Examples of X = ground-state WFs that we have used are summarized 
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in Table 3. We were unable to reproduce both binding energy and size with 
the two-parameter hyperradial function (7). In the hmit Ki Kq we get 
-Rrms = 3-0 fm. The asymptotic behaviour of the WF changes in this limit. 
In order to reproduce both binding energy and size we have used one more 
model hyperradial function 



X^^\p)^No[^—\ exp(-M (10) 



which behaves correctly both at large and small p. The strength functions, ob- 
tained with the one-parameter function ^i(^^Be) and with the two-parameter 
function of Eq. (10) ^3(^^Be), are compared in Fig. 2(b). Note that the transi- 
tion ME must be calculated numerically when (10) is used. We see that unlike 
the ^He [7] and ^^Li cases the low-energy part of the strength function, and 
thus the peak position, changes substantially when going from the one- to the 
two-parameter WF. This happens due to the fact that, at the binding energy 
and size of ■'^^Be used for fitting our bound-state WF, the non-asymptotic 
part of the hyperradial function contributes noticeably to the results. 

We also tried to vary the binding energy, the radius and the admixture of 
different components in the ground-state WF. However in all the cases we 
were unable to reproduce the narrow width and large amplitude of the exper- 
imental spectrum. A dominant K — A component in the ground-state WF 
leads to a very broad distribution of El strength with a small maximum at 
high energy (about 5Eq). When initial and final state WFs are obtained from 
model interparticlc potentials, such as the D4 potential of Ref. [5], it is also not 
possible to reproduce the shape of the spectrum. Thus inclusion of FSl does 
not solve the problem. (Note that while the D4 WF reproduces the binding 
energy and the size of ^^Be the corresponding potential leads to an unphysical 
s-wave bound state in ^^Be.) 

In Table 4 we compare the integrated strength extracted from the experimental 
data [40] with the strengths obtained with our ^3(^'^Be) WF and with the D4 
WF of Ref. [5]. When the integration is performed up to the ^^Be threshold 
the theoretical strengths prove to be substantially lower than the experimental 
one. Furthermore, an estimate of the strength integrated up to infinity can be 
made using the experimental radius, Rrms = 3. 10 ±0.15 fm [32], together with 
the El non-energy-weighted cluster sum rule 



A simple relation between -Rrms and the distance Vc between the CM and the 
core in a three-body picture is obtained by assuming (x^) = (y^) = (p^)/2 
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(see Appendix A.l of Ref. [7]). This gives 




-^rms(^) -^ri 



rms 



A 



(A) 



It is clear that the integrated strength obtained from the EMD experiment [40] 
is too large as compared with the estimates in the table. To fit the experimental 
strength one needs a larger size of ^^Be while the correlation with binding 
energy requires a smaller size. Thus we encounter one more inconsistency. 

In conclusion, we state that there are mutual disagreements between exist- 
ing experimental data on the binding energy [30,31], the radius [32], and the 
EMD energy spectrum [40] of ^^Be. In our three-body model we are unable 
to reproduce these data all at once. This result is in line with previous micro- 
scopic studies of ^"^Be [4,5,41] and in contrast with our results for other halo 
nuclei. At the same time we can confirm that a large s-wave component in 
the ^^Be ground state is required to obtain an accumulation of El strength at 
low energies. 



4 Discussion 

In Ref. [7] our analytical model was tested on the relatively well-known ^He 
nucleus. In the present paper we have studied ^^Li and ^^Be for which both ex- 
perimental data and theoretical understanding are much poorer. Concerning 
^^Li we conclude that data on El strength functions, analyzed in the frame- 
work of our model, can help to resolve the issue on the relative content of 
and components in the ground state of ^^Li. However, the exist- 

ing sets of EMD data are not consistent with each other. As our calculations 
show, the El strength function of Sackett et al. [27] seems to favor a configu- 
ration with a dominant (si/2)^ component while the data of Zinser eA, al. [26] 
and Shimoura et al. [28] are more in agreement with a mixture of and 
(pi/2)^ components with equal weights. 

Regarding ^^Be we have analyzed the available experimental data on EMD and 
found disagreements between these data, the measured size, and the binding 
energy. In the framework of our model it also proved to be impossible to repro- 
duce the narrow width of the strength function. This is in contrast with the 
results for other two-neutron halo nuclei. A theoretical study of ^^Be based on 
model intercluster potentials led to a similar disagreement [5]. This disagree- 
ment is thus not due to the neglection of FSI. The considerable weight of the 
{P1/2Y component in ^"^Be WF, suggested in [41], will also not help to reduce 
the width in our model. There is a possibility that the structure of ^^Be is more 
complicated than for other Borromean, two-neutron halo nuclei. Perhaps ^^Be 
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has more in common with ^He than with ^He in which case calculations should 
be done in the framework of a five-body model or should include excited-core 
configurations. More complete experimental data on neutron-removal cross 
sections could help to clarify whether ^^Be has a three-cluster or five-cluster 
structure. The five-body analysis of the experimental data on interaction cross 
sections, similar to that performed for ^He [42], would give a slightly smaller 
size of ^^Be than the size obtained in the framework of a three-body pic- 
ture [32]. However, in a five-body model the El strength would probably have 
an even broader distribution still not allowing the reproduction of the exper- 
imental EMD data of Labiche et al. [40]. Furthermore, methods to measure 
masses of short-lived isotopes have improved significantly since the experi- 
ments [30,31] were performed, and new experimental results on the ^^Be mass 
would be very important. 
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Fig. 1. A comparison of our strength function with experimental data for ^^Li (cir- 
cles [26], triangles [27], squares [28]), both in arbitrary units (a) and in absolute scale 
(b). The model WFs are: dotted - *i(^^Li), dashed - *2(^^Li), solid - *3(^^Li), 
dash-dotted - *4(^^Li), see Table 1. 
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Fig. 2. The positions of the peak obtained with two different binding energies of ^^Be 
(a). The curves correspond to the one-parameter function (6) with Eq = 1.34 MeV 
(dash-dotted) and Eq = 0.90 MeV (dotted). In (b) it is shown how the distribution 
changes when trying to fit both binding energy and size. Here the solid curve corre- 
sponds to the two-parameter *3(-'^^Be) WF, Eq. (10), which reproduces both these 
quantities, while the dash-dotted curve is the same as in (a). In both figures a com- 
parison is made with "experimental data" obtained by dividing out the spectrum 
of virtual photons from the energy spectrum of Ref. [40]. 
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Table 1 

The model WFs for ^^Li. To reproduce binding energy and rms radius we use 
kq = 0.1193 fm~^ and ki = 0.3173 fm~^. The quantity Tc is the distance between the 
core and the CM. The approximate amplitudes of and configurations 

of valence neutrons are shown in the last two columns. 



WF 


x(p) 


«00 


020 


021 


i?i-ms 


rc 








Eq. 








(fm) 


(fm) 






*i(^^Li) 


(6) 


1 








2.76 


0.60 


1 





*2(''Li) 


(7) 


1 








3.55 


0.95 


1 





*3(''Li) 


(7) 








3.55 


0.80 


Vol 


Vol 


*4(''Li) 


(7) 




VV6 




3.55 


1.08 


Vol 


-Vol 



Table 2 

Comparison of hyperradii and binding energies of three Borromean halo nuclei. 



Nucleus £^0 (MeV) p^ms (fm) 



iiLi 


0.295 


9.5 


^He 


0.97 


5.4 


i^Be 


1.34 


7.4 



Table 3 

The model WFs for ^^Bc. The binding energy of = 1.34 MeV corresponds to 
Kq = 0.2543 fm-i while Eq = 0.90 MeV corresponds to = 0.2084 fm^^. The 
parameter po entering the ^'3(^'*Be) WF, Eq. (10), equals 5.42 fm. The quantity 
Tc is the distance between the core and the CM. Only the K = term in the HH 
expansion has been retained. 



WF 


x{p) 


Eq 








Eq. 


(MeV) 


(fm) 


(fm) 


*i(i4Be) 


(6) 


1.34 


2.51 


0.21 


*2(''Be) 


(6) 


0.90 


2.56 


0.26 


*3('^Be) 


(10) 


1.34 


3.10 


0.57 
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Table 4 

Experimental and theoretical integrated El strengths. The first row represents 
strengths integrated up to 3.2 MeV (^^Bc threshold). In the second row total 
strengths are shown. The total strength in the second column is calculated from 
the value -Rrms = 3.10 it 0.15 fm [32] using the non-energy-weighted cluster sum 
rule. To this aim, the Tc value of 0.57 it 0.07 fm was obtained from the above value 
of -Rrms assuming a simple three-body structure. 





Extracted from [40] Obtained from -Rrms [32] 


*3(^'Be) 


D4 [5] 


5e<3.2(E1) (e2fm2) 


1.40 ± 0.40 


0.70 


0.96 


5(E1) {eHxn^) 


1.23 ± 0.33 


1.23 


1.56 
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